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ABSTRACT 

We compare the theoretical predictions of the transition layer model with some observational features of 
quasi-periodic oscillations (QPOs) in neutron star X-ray binaries. We found that the correlation between 
horizontal branch oscillation (HBO) frequencies and kilohcrtz (kHz) QPO frequencies, the difference 
between the low-frequency QPOs in atoll sources and HBOs in Z sources, and the correlation between 
the frequencies of low- frequency QPOs and break frequencies can be well explained by the transition layer 
model, provided the neutron star mass is around 1.4 solar mass and the angle between magnetosphere 
equator and accretion disk plane is around 6 degree. The observed decrease of peak separation between 
two kHz QPO frequencies with the increase of kHz QPO frequencies and the increase of QPO frequencies 
with the increase of inferred mass accretion rate are also consistent with the theoretical predictions of 
transition layer model. In addition, we derive a simple equation that can be adopted to estimate the angle 
(5) between magnetosphere equator and accretion disk plane by use of the simultaneously observed QPO 
frequency data. We estimate these angles, in the range of 4 to 8 degrees, for five Z sources and two atoll 
sources. The nearly constant S value for each source, derived from the different sets of simultaneously 
observed QPO frequency data, provides a strong test of the theoretical model. Finally, we suggest that 
the similar transition layer oscillations may be also responsible for the observed QPOs in accretion- 
powered millisecond X-ray pulsar and Galactic black hole candidates. 

Subject headings: accretion, accretion disks — stars: neutron — stars: binaries: general — X-Rays: 
stars 



1. INTRODUCTION 

Quasi-periodic Oscillations (QPOs) have been found to 
be a very common feature of accreting systems around 
compact objects. In neutron star low-mass X-ray bina- 
ries (LMXBs) with low magnetic field, at least four types 
of QPOs, namely the 15-60 Hz horizontal branch oscilla- 
tions (HBOs; van der Khs et al. 1985), the 6-20 Hz nor- 
mal branch oscillations (NBOs; Middlcdich & Pricdhorsky 
1986) and the 300-1200 Hz lower and upper kilohertz (kHz) 
QPOs (van der Klis et al. 1996), have been discovered by 
several X-ray satellites including EXOSAT, GINGA, and 
the Rossi X-Ray Timing Explorer (RXTE). In spite of the 
different evolution tracks in the X-ray color-color diagram 
(Hasinger & van der Klis 1989), the low- frequency QPOs 
in atoll sources seem to have the similar feature to those of 
the HBOs in Z sources (Strohmayer et al. 1996, Wijnands 
& van der Klis 1997). With the advantage of RXTE, 20 
sources have now shown kHz QPOs and 18 of them have 
shown two simultaneous kHz peaks (see van der Klis 2000 
for a review). The peak separations of two kHz QPOs 
have been observed to decrease considerably in five sources 
when the kHz QPO frequencies increase (van der Klis et al. 
1997; Mendez et al. 1998a,b; Ford et al. 1998, Mendez & 
van der Klis 1999, Markwardt et al. 1999). Other sources 
may also show similar variations (Psaltis et al. 1998). It 
has been found that the HBO frequencies of Z sources and 
the frequencies of low-frequency QPOs in atoll sources are 
well correlated with the kHz QPO frequencies (Psaltis et 
al. 1999a). Similar correlation of low and high frequency 
QPOs might also exist for some other neutron star systems 



and Galactic black hole candidates. In addition, tight cor- 
relations of the HBO frequencies of Z sources and of the 
frequencies of low-frequency QPOs in atoll sources and 
black hole candidates with the break frequencies shown in 
their power density spectra have also been found recently 
(Wijnands & van der Klis 1999). These two correlations 
strongly imply that similar mechanisms may be responsi- 
ble for the break frequency, the low-frequency QPOs and 
the high-frequency QPOs in both neutron star and black 
hole X-ray binaries. 

There are several QPO models that are currently in de- 
bate. Due to the difficulty in producing kHz QPOs, the 
previous magnetospheric beat-frequency model proposed 
by Alpar & Shaham (1985) has been modified as sonic- 
point beat-frequency model (Miller, Lamb & Psaltis 1998). 
In this model the upper kHz QPOs arise from the clumps 
moving with the Keplerian frequency at the sonic point 
that is near the innermost stable circular orbit of neutron 
star. The lower kHz QPO frequency is thus the beat fre- 
quency between the upper kHz QPO frequency and the 
neutron star spin frequency. However, this model, if with- 
out further modifications, predicts a constant peak separa- 
tion between two kHz QPO frequencies, which is contrary 
to the observations. Moreover, recent investigation on the 
boundary layer accretion onto neutron star indicated that 
the accretion flow is always subsonic and there is probably 
no such a sonic point at all near the neutron star surface 
(Popham & Sunyaev 2001). This leads to further doubt 
about the sonic-point beat-frequency model. Another at- 
tractive model is the relativistic precession model(RPM) 
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proposed by Stella & Vietri (1998, 1999), who identified 
the upper kHz QPO frequency with that of an Kcplerian 
orbit in the disk and the lower kHz QPO frequency and 
HBO frequency with, respectively, the periastron preces- 
sion frequency and twice of the nodal precession frequency 
of this orbit. This model can qualitatively explain the ob- 
served correlations between these three QPO frequencies 
(Stella, Vietri & Morskink 1999), but the precise match 
requires some fine tuning of additional free parameters 
such as the orbital eccentricity and the periastron dis- 
tance. Moreover, in order to explain the HBOs and its 
correlation with kHz QPOs, this model requires that the 
I/M value (where I is the momentum of inertial and M is 
the mass of neutron star) be 2 times larger than that pre- 
dicted by realistic neutron star equation of state (Psaltis 
et al. 1999b). In addition, RPM model requires the neu- 
tron star mass in the range of 1.9-2.2M0 in order to ex- 
plain the correlation between low and high QPO frequen- 
cies, which is significantly larger than the measured value, 
l.lAfQ<M<1.6MQ, for neutron stars in radio pidsar bi- 
naries (Thorsett & Chakrabarty 1999; Finn 1994). The 
problem to make a neutron star with mass larger than 
2.2Mq arises because it would require accretion of mate- 
rial of at least 0.8Mq, which means that these X-ray bi- 
naries would have to be rather old and would have to have 
spun up rapidly for 90% of their lifetime (Lai, Lovelace 
& Wasserman 1999). The same problem probably exists 
for the new model proposed by Psaltis & Norman (2000), 
who derived the similar characteristic QPO frequencies as 
in RPM in the inner accretion disk. 

Another alternative model for QPOs in neutron star 
X-ray binaries is the disk transition layer model (hear- 
after TLM) proposed recently by Titarchuk & Oshcrovich 
(1999), Osherovich & Titarchuk (1999a,b) and Titarchuk, 
Oshervich & Kuznetsov (1999). The geometry of the tran- 
sition layer model has been clearly shown in Figure 1 of 
Titarchuk et al. (1999). In this "two-oscillator" model, the 
lower kHz QPO frequency is identified with the Kcplerian 
frequency in a viscous transition layer between Kcplerian 
disk and neutron star surface. Assuming the magnetic field 
of neutron star is low (in the range of 10^ to 10^ Gauss), 
the Alfvcn radius, depending mainly on the magnetic field 
and accretion rate, is at most of several radii of neutron 
star. The size of the magnetosphere, as well as that of the 
transition layer, is thus very small. Viscous oscillations 
in this transition layer produce the observed noise break 
and a low frequency QPO. In addition, the adjustment 
of the Kcplerian disk to the sub-Keplerian layer may cre- 
ate conditions favorable for the formation of a hot blob in 
the transition layer. This blob, when thrown out into the 
rotating magnetosphere from the transition layer, partic- 
ipates in the radial oscillations with Kcplerian frequency. 
Under the infiuence of the Coriolis force, such a blob, as- 
sumed to be a Kcplerian oscillator, oscillates both radially 
and perpendicular to the disk. This produces two harmon- 
ics of another low-frequency QPO (HBO in Z sources) and 
the upper kHz QPO. These six different frequencies have 
been identified in two atoll sources (4U 1728-24 and 4U 
1702-42) and Z source Sco X-1 ( Titarchuk et al. 1999; 
Osherovich & Titarchuk 1999a, b). The observed corre- 
lations of QPO frequencies for these sources seem to be 
consistent with the theory of TLM. 

In this paper, we derive more theoretical predictions 



for TLM and compare them with the available observa- 
tional data of QPOs in neutron star X-ray binaries. The 
consistency between model predictions and observational 
data without assuming larger neutron star mass and a stiff 
equation of state suggests that TLM is perhaps a more 
competitive model than others. We also suggest that a 
similar model may be also applied to the cases of millisec- 
ond X-ray pulsars and Galactic black hole candidates. 

2. CHARACTERISTIC QPO FREQUENCIES IN THE 
TRANSITION LAYER MODEL 

As indicated by Titarchuk & Osherovich (1999) and Os- 
herovich & Titarchuk (1999b), there are six characteristic 
QPO frequencies in the TLM. First, the fundamental fre- 
quency is assumed to be the Kcplerian frequency, namely, 

where G is the gravitational constant, M is the mass of 
neutron star and R is the radius of an orbit in the transi- 
tion layer. The linear Kcplerian oscdllator with frequency 
f2fc/27r in the frame of reference rotating with rotational 
frequency 0/27r (not perpendicular to the disk plane) is 
known to have an exact solution describing two branches 
of oscillations (see §39 in Landau & lifshitz 1960). Assum- 
ing fl is constant, the dispersion relation of the frequency 
of oscillation uj can be derived as: 

uj^ - {fll + 402)a;2 + m'^nlsin'^S = 0, (2) 

where 5 is the angle between fi and the normal to the Kc- 
plerian oscillation. If S is assumed to be small, one can 
obtain simple expressions of two characteristic oscillation 
frequencies from the dispersion relation. The radial eigen- 
mode has a frequency 

= sj^ + i^?^ (3) 

where Q. is the angular velocity of the rotating magneto- 
sphere in the case of neutron star X-ray binaries. If the 
magnetosphere corotates with the neutron star, we have 
= f2o where f^o is the rotating angular velocity of neu- 
tron star. The frequency Vh is the analog of the hybrid 
frequency in plasma physics (Akhiezer et al. 1975; Benson 
1977). On the other hand, the frequency of the vertical 
eigenmode is 

vl = ——sinS, (4) 

These three characteristic frequencies, j/fe, i^h and j/^, are 
assumed to account for the frequencies of lower kHz QPOs, 
upper kHz QPOs and HBOs in Z sources, respectively. Be- 
cause the second harmonics of HBO are often seen in Z 
sources, they are interpreted as 2j/l, counting as the other 
characteristic frequency. 

From equations (3) and (4), we can derive several very 
useful relations among these characteristic frequencies. 
First; the correlation between the HBO frequency and two 
kHz QPO frequencies can be expressed as: 



i^L = i'k\n-{ — Ysin5. (5) 
V 
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Second, the difference between two kHz QPO frequencies 
can be related to as 



Az/ = z/ft - i/fe = vk{\j 1 + { — f 



!)• 



(6) 



Tliird, if tlie Keplerian oscillation is aligned with the accre- 
tion disk plane, the angle between magnetosphere equator 
and disk plane can be determined by 



= arcstn{ ====). 



(7) 



Because three frequencies in the right of equation (7) are 
all observational quantities, this equation perhaps suggests 
a useful method to determine 6 for neutron star X-ray bi- 
naries. 

The other two characteristic frequencies in the TLM 
are related to the two typical timescalcs in the transition 
layer. One is the radial drift timescale of matter through 
the layer bounded between disk and neutron star, namely, 
tr ~ {Rout — Ro)/vr, where Rq is the radius of neutron 
star (taken to be GGM/c^ in this paper) and is the ra- 
dial drift velocity. This timescale gives the characteristic 
viscous frequency as 



71/ 



{Rout — Ro)R 
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R rout - 1 



(8) 



where rout is defined as rout = Rout/Ro, and 7 is the 
Reynolds number given by 7 = Ad / AiTpvH = VrRjv 
where M is the mass accretion rate and p is the mass 
density. We adopt the standard viscosity prescription 
V = aCsH = aVikH"^ : where Cs is the local soimd speed 
and H is the scale height (Shakura & Sunyaov 1973). The 
relation between 7 and Tout can be obtained by considering 
the radial transport of angular momentum and the bound- 
ary conditions at i?o and Rout (Titarchuk & Osherovich 
1999). Assuming that at i?o, = fio and at Roui, = f^fc 
and dfl/dr = dflk/dr and defining A — i}ko/^o (where 
rifco is the Keplerian angular velocity at Rq), rout can be 
determined by following equation: 
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j)Ari-,^ + {^-2)r;Z + ^A = 0. 



(9) 



Another characteristic timescale in the transition layer 
is the diffusion time scale given by ^d?// ^ {Rout — 
RqY /{IdVr), where Id is the diffusion length scale in the 
transition layer. It gives the break frequency as 



^break 



Id 



Id 



{Rout — Ro) Ro Tout — 1 



(10) 



This expression of break frequency is similar as that given 
in equation (29) of Psaltis & Norman (2000) if Id is re- 
placed by R. Similar to the argument given by Psaltis 
& Norman (2000), the reason that Vhreak observed as a 
frequency break may be that the transition layer acts as 
a filter band so that the response is constant if the oscil- 
lation frequency is lower than the inverse diffusion time 
scale and the response decreases above it. Titarchuk & 
Osherovich (1999) described Id, as the mean free path of a 



particle, however, this introduced another uncertain quan- 
tity in the TLM theory. In this paper, using a similar 
approach adopted in the standard accretion disk model 
(Shakura & Sunyaev 1973), we assume that Id is smaller 
than the vertical and radial size of the transition layer and 
take Id/Ro as a constant less than unity for simplicity (in 
§3.3 we will show this assumption is consistent with the 
observed data of break frequency). The equations (8), (9) 
and (10) can be used to calculate the correlations between 
other low-frequency QPO (besides HBO) frequency, break 
frequency and Keplerian frequency. The original descrip- 
tions of viscous and break frequencies shown in equations 
(12) and (13) in Titarchuk & Osherovich (1999) have been 
improved here (see equations (8) and (10)) by considering 
a more explicit expression of the radial drift timescale. 

3. COMPARISONS OF MODEL PREDICTIONS WITH THE 
OBSERVATIONAL DATA 

There are a lot of observations made by several X-ray 
satellites on QPOs in neutron star X-ray binaries. But 
only until very recently the correlations among the fre- 
quencies of low and high-frequency QPOs, the variation 
of kHz QPO frequencies and the difference between HBOs 
and other low-frequency QPOs have been systematically 
studied with plenty of QPO data. In this section we try 
to compare the predictions from the TLM, as indicated 
by equations in the section above, with the observational 
results. 

3.1. Correlation between HBOs and kilohertz QPOs in Z 

sources 

Recently, Psaltis et al. (1999a) have compiled the 
low and high-frequency QPO data of Z sources and 
atoll sources, as well as black hole candidates. They 
found a tight correlation between HBO frequency and 
the frequency of lower kHz QPO. The data points for 
the Z sources arc consistent with the empirical relation 
fHBO ~ (42 ± 3Hz)(i/i/(500Hz))0-95±o i6when the fre- 
quency of lower kHz QPO Pi less than 550 Hz. When ui 
larger than 500 Hz, vrbo increases slowly with vi. This 
feature seems to be difficult to fit by the RPM. 

In Figurel, we plot vhbo and vi data for five Z sources 
and compare it with the predictions from both RPM 
and TLM. The fit parameters for RPM (M = 1.95M0, 
a/M = 0.22) are taken from Stella et al. (1999) and a is 
the spin parameter of the neutron star. It is evident that 
RPM predicts a steeper correlation than the observed one 
and can not explain the flatness of such a correlation when 
U\ > 500 Hz. Any fine tuning of the fit parameters M and 
a/M can not provide a better fit. However, when we plot 
the theoretical curves of TLM according to the equation 
(5) by choosing quite reasonable parameters Vh — ^k = 300 
Hz and 5 around 6°, the match of model prediction and 
observation data is quite well. In addition, we also plot 
the observed correlation between vhbo and upper kHz 
frequency 1^2 for five Z sources in Figure 2. In this case, 
the RPM can not fit the observation data if we still choose 
the parameter M = 1.95M0, a/M = 0.22. Even when we 
choose the fit parameter M = 1.95Af0, a/M = 0.30, the fit 
is still not satisfactory. Therefore, our comparison shows 
that RPM, at least in its current version, is difficult to 
fit the correlations between simultaneously observed HBO 
frequencies and kHz QPO frequencies. On the contrary, 
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Fig. 1. — The correlation between HBO frequencies and lower kHz QPO frequencies of Z sources. The solid, short-dashed and long-dashed 
lines represent TLM predictions with parameters \v = 300 Hz and <5 = 5.5°, 6° and 6.5°, respectively. The dotted line represents the 
prediction of RPM with M = 1.95Mq and a/M = 0.22. The error bars, if not shown, are comparable to the size of the symbols. 



from Figure 1 and Figure 2 we can see clearly that TLM 

predicts rather flatter slopes than RPM. With the same 
parameters, TLM can provide good fits to both the I'm bo 
— ui and vhbo ~ ^2 correlations. 

For some sources (especially GX17+2 and GX 5-1) in 
Figurel, it seems that the slope is steeper than that pre- 
dicted by TLM and perhaps matches better with the RPM 
prediction. However, we must note that the TLM predic- 
tions plotted in figures 1 and 2 were based on a very simple 
assumption {vh—fk = 300Hz). In fact, the peak separation 
between two kHz QPO frequencies, Vh — ^k, is certainly not 
constant for individual sources (see §3.4). If we take the 
observed values of Vh and directly and use equation (5) 
to predict vl, the agreement with the observations will be 
more close for individual sources plotted in figures 1 and 2. 
On the other hand, with reasonable parameters, we noted 
that RPM is unable to fit both the observed correlations 
shown in figures 1 and 2 consistently. Even if the RPM 
prediction could better match the slopes of 0X17+2 and 
GX 5-1 in Figure 1, it fails to match the observed slopes 
for these sources in Figure 2 if the same parameters were 
used. 

3.2. Difference between low frequency QPOs in atoll 
sources and HBOs in Z sources 



Although there are some reports on the similar QPOs 

in atoll sources to HBOs in Z sources (Strohmayer et al. 
1996, Wijnands & van der Klis 1997), it can be clearly seen 
from Figure 2 in Psaltis et al. (1999a) that the frequencies 
of these low-frequency QPOs (10-50 Hz) of atoll sources 
have steeper dependence on the lower kHz QPO frequency 
than those of HBOs in Z sources. In Figure 3, we plot 
this correlation and compare it with the empirical relation 
found for HBOs in Z sources. Obviously the difference be- 
tween these low-frequency QPOs and HBOs is large. It 
is therefore quite possible that these low-frequency QPOs 
may have different origin from HBOs in Z sources. Actu- 
ally, the similar low-frequency QPOs (usually seen as extra 
noise components) have been found simultaneously with 
HBOs in Z source Sco X-1 and GX 17-1-2 (van der Klis et 
al. 1997; Titarchuk et al. 1999; Wijnands & van der Klis 
1999), while the similar HBO-likc frequencies have been 
identified in two atoll sources 4U 1728-34 and 4U 1702-42 
(Titarchuk et al. 1999; Osherovich & Titarchuk 1999b). 

In Figure 3, it can be clearly seen that these low- 
frequency QPOs in Z source Sco X-1 follow the same cor- 
relation with vx as other atoll sources. This probably sug- 
gests that the low- frequency QPOs in both Z and atoll 
sources have the same physics origin and they are obvi- 
ously different from HBOs. In TLM, the frequency of this 
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Fig. 2. — The correlation between HBO frequencies and upper kHz QPO frequencies of Z sources. The hnes have the same meanings as in 
Figurel. The dot-dashed Une represents the prediction of RPM with M = l.gSM© and a/M = 0.30. 



kind of QPOs is described as the viscous frequency in the 
transition layer, therefore the mechanism to produce the 
low-frcqucncy QPOs is clearly different from HBOs. Using 
equations (8) and (9) and assuming a{H/Rf = 7.0 x 10"'', 
M around IAMq and Tout in the range of 1.2 to 1.8 
(namely 7 in the range of 6 to 20), we plot the model pre- 
diction of TLM in Figure 3 (the result is insensitive to the 
spin frequency of neutron star, z/q = 300 Hz is assumed to 
make the plot). We can see the close agreement between 
theory and observations. For KS 1731-26 and 4U 1735-44, 
the data points are a little away from others. They are 
likely associated with the NBOs of Sco X-1 (Psaltis et al. 
1999a), though they can be still fit by assuming a smaller 
value of a{H/Kf. However, we noted that for 4U 1735-44, 
recent observations indicated the existence of a 67Hz QPO 
together with a possible 900Hz lower kHz QPO (Wijnands 
et al. 1998c). If we plot this point in Figure 3, it will also 
locate in the range of TLM prediction. 

In Figure 3, only the data of February 16, 1996 for 4U 
1728-34 are plotted because in these observations the low 
kHz QPOs can be clearly distinguished from the simulta- 
neously observed upper kHz QPOs (Ford & van der Klis 
1998). The agreement of the observation data with our 
predictions assuming M around 1.4Mq suggests that the 
compact object in 4U 1728-34 may not be a strange star 
as argued by Li et al. (1999). We note that their conclu- 



sion is based on the fit result, ak = 1.03, by Titarchuk & 
Osherovich (1999), who used less accurate expression for 
viscous frequency (sec §2) and ambiguous data for lower 
kHz QPO frequencies of 4U 1728-34. 

3.3. Correlation between frequencies of low-frequency 
QPOs and break frequencies in Z sources and atoll 
sources 

A recent detailed study on the broadband power density 
spectra of X-ray binaries indicates that the frequency (1- 
60 Hz) of low-frequency QPOs (including HBOs) have a 
strong correlation with the break frequency [vbreak] 0.1-30 
Hz) of both Z sources and atoll sources, possibly as well 
as black hole candidates (Wijnands & van der Klis 1999). 
However, the HBO frequencies of Z sources follow a differ- 
ent relationship with Vbreak from the low-frequency QPOs 
of other sources. Such a difference has not been satisfacto- 
rily explained so far (see, however, Titarchuk et al. (1999) 
for a possible explanation but note that they mistakenly 
plotted the HBO frequencies of Z sources as the viscous 
frequencies in their Figure 4). 

In TLM, because HBO and viscous frequencies come 
from different mechanisms, it is natural that they follow 
different correlations with the break frequencies. Accord- 
ing to equation (5) and equations (8)-(10), we can derive 
the predicted correlation of HBO frequencies, viscous fre- 
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Fig. 3. — The correlation between the frequencies of low-frequency QPOs and lower kHz QPO frequencies. The solid, short dashed and 
long dashed lines represent the TLM predictions with a{H/K)'^ = 7.0 X 10""' and M =1.2Mq, IAMq and LOA/q, respectively. The dotted 
line represent the prediction with a(H/R)^ = 3.2 X 10~^ and M =1.4Mq. The dot-dashed line represents the empirical relation between 
HBO frequencies and lower kHz QPO frequencies for Z sources. 



quencies with break frequencies. In Figure 4, we plot these 
predictions and compare them with the observation data. 

For atoll sources, our prediction closely agrees with the 
relation, Vbreak = 0.04j^]^|,'^, which was derived from the 
observed correlations vlp - VkHz and Vbreak i^kHz for 
atoll source 4U 1728-34 (Ford & van der Klis 1998). The 
viscous frequencies of other atoll sources and Z source Sco 
X-1 also follow nearly the same v^p fbreak correlation. 
We adopted Id/Ro = 0.3 in equation (9) and the same 
parameters {a{H/Rf' = 7.0 x 10~^, M=1.4M0) as above 
to derive this correlation. The value, Id/Ra = 0.3, is ap- 
preciated by the observed data and also consistent with 
our assumption that the diffusion length scale is smaller 
than the vertical and radial size of the transition layer. 
Note that = O.SRq < H means H/R may be larger 
than 0.3. This is quite possible if the magnetic field is 
low and the radial size of transition layer is small. If we 
take a{H/R)^ = 7.0 x 10~^, it requires that a should be 
less than 0.01. Such a lower value of a is still well within 
the range of viscosity parameter discussed in the accretion 
disk model with boundary layer. Assuming Vh — Vk = 300 
Hz and 5 around 6°, we derive the correlation between 
HBO frequencies and break frequencies. Comparing with 
the observation data of Z sources, we note that the agree- 



ment is well. There are several data points for atoll source 
4U 1728-34 locating at the branch for HBO in Z sources, 

which may indicate that these data are actually frequen- 
cies of HBO-typc QPOs in atoll sources (see §3.2). 

3.4. Variation of peak separation of two kilohertz QPOs 

The observed variation of difference of two kHz QPO 
frequencies in at least one Z source, Sco X-1 (van der Klis 
et al. 1997) and four atoU sources, 4U 1728-34, 4U1608-52, 
4U 1735-44 and 4U 1702-43 (Mendez et al. 1998a,b; Ford 
et al. 1998; Mendez & van der Klis 1999; Markwardt et 
al. 1999), seem to exclude any QPO model that predicts 
the constant peak separation. The magnetospheric beat- 
frequency model, as well as the sonic-point beat-frequency 
model, probably can not be responsible for the kHz QPOs, 
unless further modifications are made. 

Recent precise measurements of kHz QPO frequencies 
show that the peak separation Au = 1/2 — 1^1 decreases 
considerably when kHz QPO frequencies increase (van der 
Klis 2000). This feature can not be well fitted by the RPM, 
which predicts a much steeper decrease of Ai> when 1/2 in- 
creases to more than 1000 Hz and a significant decrease of 
Al/ when 1^2 decreases to less than 700 Hz (Stella & Vietri 
1999). The latter effect has not been observed yet. A pre- 
cise match between RPM and observations requires some 
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Flc;. 4. — The correlation between the frequencies of low-frequency QPOs and break frequencies. The short dashed, long dashed and 
dot-dashed lines represent the TLM predictions with At/ = 300Hz and <5 = 4.5°, 6°, 7.5°, respectively. The solid line represent the prediction 
with a{H/Rf = 7.0 x lO"* and M =IAMq. I^/Rq = 0.3 was assumed to derive the theoretical break frequency. The dotted line represents 
the empirical relation inbreak = 0.04i/|^|;^ derived for 4U 1728-34. Other atoll sources are plotted as open triangles, diamonds and squares. 



additional assumptions such as the dependence of orbital 
eccentricity on the orbital frequency. In TLM, however, 
the decrease of Au as the increase of kHz QPO frequency 
is a natural result even if we simply take fl = constant 
(see equation (6)). In Figure 5, we show the Av — v\ rela- 
tion for all 5 sources with observed Av variations. When 
we plot the TLM prediction assuming the constant 0/27r 
(340Hz, 365Hz and 380Hz) in equation (6), we see that 
even at this simple assumption the theoretical predictions 
are closely agreement with most observation data of Sco 
X-1, 4U 1708-52, 4U 1735-44 and 4U 1702-43. However, 
it is clear that this simple assumption can not explain 
the data of 4U 1728-34. In fact, the matter in magne- 
tosphere probably experiences differential rotation and 
may not be a constant. If we assume that the angular 

velocity could be expressed ^ = vq — {civ^^^ — C2i'^)'^ 
(Osherovich & Titarchuk 1999a) where vq, ci and C2 are 
constants, and take uq = 405 Hz, ci = 0.26Hz~^/^ and 
C2 — 3.4 X 10~'''Hz~^/^, we can sec the prediction of peak 
variation of twin kHz QPO frequencies is close to the ob- 
servation data for 4U 1728-34. In summary, it seems that 
the variation of peak separation of two kHz QPOs is also 
consistent with the theoretical prediction of TLM. 



4. OBSERVATIONAL DETERMINATION OF MISALIGNMENT 
OF NEUTRON STAR MAGNETOSPHERE AND 
ACCRETION DISK 

There are only few methods that can be used to es- 
timate the magnetic inclination angle between magnetic 
axis and neutron star rotation axis for radio pulsars (Lyne 
& Manchester 1988). It may be more difficult to measure 
these angles for X-ray pulsars and low-mass X-ray bina- 
ries. However, according to TLM (see equation (7)), if we 
have simultaneous observation data for frequencies of both 
kHz QPOs and HBOs, wc may estimate the angle between 
magnetopherc equator and disk plane for neutron star X- 
ray binaries. We note that equation (7) is derived under 
the assumption that 5 is small. In fact, a more accurate 
equation, which is also applicable to arbitrary values of 
(5, can be directly obtained from the dispersion relation of 
Keplerian oscillator (see equation (2)), namely 



5 = arcsin{ 



V2VHB0 



v\\fv. 



(11) 



'HBO 



where vx^ vi and vhbo are the observed lower, upper kHz 
QPO and HBO frequencies respectively, and we assume 
that they correspond to the Keplerian frequency and two 
characteristic oscillation frequencies derived from the dis- 



8 




J I I I I I I I I I I I I I I I I I I I I I I I I L 

500 600 700 800 900 

Lower KHz QPO frequency (Hz) 



Flc;. 5. — The relation between the peak separations of two kHz QPO frequencies and the lower kHz QPO frequencies. The solid, dotted 
and short dashed lines represent the TLM predictions using equation (6) with VIq/2tt =340 Hz, 365 Hz and 380 Hz, respectively. The long 
dashed line represent the prediction using equations (3) and (6) with parameters vq = 405 Hz, ci = 0.26 Hz"^/*" and C2 = 3.4 x 10~^ Hz~^/-^. 



persion relation (see §2). In the limit uhbo « t^2, equa- 
tion (11) is identical to equation (7). Note that equation 
(11) can be also used as a test of TLM because the derived 
values of S using dilTerent set of simultaneously observed 
QPO frequency data for each source should be nearly the 
same (Titarchuk & Osherovich 2000). 

Up to know, simultaneous frequency data of both kHz 
QPOs and HBOs have been obtained with RXTE for five 
Z sources. In two atoll sources with simultaneous data of 
two KHz QPO frequencies, the HBO-like frequencies have 
been identified recently (Titarchuk et al. 1999; Osherovich 
& Titarchuk 1999b). Therefore, the misalignment of mag- 
netosphere and accretion disk can be estimated for these 
sources according to equation (11). In table 1, we list the 
observation data and the derived average values of S for 
them. It is clear that all of them have smaller 6 values and 
are therefore consistent with the assumption in §2. Espe- 
cially for Z sources, their S values are in a very narrow 
range from 5. "4 to 6. "4. Moreover, the standard deriva- 
tions of derived S values are small, which means nearly 
the same 6 value is obtained for each source using the si- 
multaneously observed QPO frequency data. This is well 
agreement with the prediction of TLM. In addition, we 
noted that the method we used to determine 5 is better 
than that used by Osherovich & Titarchuk (1999b) and 



Titarchuk et al. (1999), who derived S by fitting fl calcu- 
lated with the observed kHz QPO frequencies (see equa- 
tion (2)). It is obviously not necessary to do so if we have 
simultaneously observed frequency data of HBOs and kHz 
QPOs. Instead, it is more straightforward and accurate to 
derive S by using equation (11). However, their estimated 
values, S = 3.°9 ± 0.°2 for 4U 1702-42, S = 8.°3 ± l.°0 for 
4U 1728-34 and S = 5.°5 ± 0.°5 for Sco X-1, are very close 
to our estimations. For other four Z sources in Table 1, 
it is the first time to give the estimations of their angles 
between magnetosphere equator and disk plane. 

5. DISCUSSION 

The comparisons of the model predictions of transition 
layer model with the observed data seem to suggest that 
TLM may provide better explanations on them than the 
current version of other QPO models. In contrary to the 
inability of beat-frequency model to explain the variation 
of peak separation of two kHz QPO frequencies and the 
difficulty of relativistic precession model to explain both 
the correlations of HBO frequencies with lower and upper 
kHz QPO frequencies, TLM is so far the only model that 
can self-consistently explain these observed correlations. 
In addition, TLM suggests some physics mechanisms for 
producing low-frequency QPOs and break frequencies in 
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Table 1 

Angle between Magnetosphere equator and accretion disk plane determined by simultaneously observed QPO 

frequency data 



Name 


Type 


^2 




vhbo 


Ref. 


<5> 






(Hz) 


(Hz) 


(Hz) 




n 


Cyg X-2 


Z 


856 


532 


46 


1,2 


6.34+0.14 


GX 17+2 


z 


640-1090 


480-780 


25-61 


1,3 


6.30+0.37 


Sco X-1 


z 


870-1080 


550-860 


41-48 


1.4 


5.44+0.22 


GX 5-1 


z 


500-890 


210-660 


17-51 


1,5 


5.80+0.62 


GX 340+0 


z 


560-820 


250-500 


22-45 


1,6 


6.40+0.45 


4U 1728-34 


Atoll 


930-1130 


600-790 


52-90 


1,7,8 


8.01+1.06 


4U 1702-42 


Atoll 


1000-1080 


660-770 


33-40 


9,10 


3.96+0.23 



References: (1). Psaltis et al. (1999a); (2). Wijnands et al. 1998a; (3). Wijnands et al. 1997; (4) van der Klis et al. 
1997; (5). Wijnands ct al. 1998b: (6). Jonkcr ct al. 1998: (7). Strohmayer et al 1996; (8). Ford & van der Klis 1998; (9). 
Markwardt et al. (1999); (10). Osherovich & Titarchuk (1999b) 



X-ray binaries, and can explain the observed correlation 
between the frequencies of low-frequency QPOs and the 
break frequencies. On the contrary, no clear mechanisms 
for these low frequency features have been discussed yet 
in most of other QPO models. 

Another clear difference between TLM and other mod- 
els is that the lower kHz QPO frequency, rather than the 
upper kHz QPO frequency, is described as Keplerian fre- 
quency. Therefore, it is not necessary in TLM to involve 
larger neutron star mass and stiff equation of state as re- 
quired in other models such as sonic-point beat-frequency 
model (Miller et al. 1998) and rolativistic precession model 
(Stella & Vietri 1998, 1999). The compact star in TLM can 
be still a normal neutron star with mass around 1.4M0 and 
radius around 12km. Indeed, the observed correlation of 
QPO frequencies can be well explained by taking these pa- 
rameters in TLM. A larger neutron star mass (> ^Mq) is 
not appreciated since it will lead to steeper correlation of 
low-frequency QPO frequencies with lower kHz QPO fre- 
quencies than that observed for atoll sources (see Figure 
2). 

The variations of the size of magnetosphere and the size 
of transition layer due to the change of mass accretion 
rate may account for the the observed variations of QPO 
frequencies in X-ray binaries. Because the size of mag- 
netosphere can be described approximately by the Alven 
radius Ra, which is inverse proportional to M"^/'', the 
magnetosphere will shrink if accretion rate M increases. 
This may lead to the increases of the kHz QPO frequen- 
cies {vh, yk) and the HBO frequency {vl)- The size of the 
transition layer can be described by Rout — Ro , which is 
directly related to the Reynolds number of accretion flow 
(sec equation (9)). The increase of accretion rate will cause 
the increase of Reynolds number, and therefore lead to the 
decrease of Rout—Ro- According to equations (8) and (10), 
both the viscous frequency and the break frequency will 
increase as the increase of accretion rate. Indeed, many 
observations have shown the increases of kHz QPO fre- 
quencies, HBO frequency and break frequency when the 
sources move along the tracks in the color-color diagram, 
equivalently, in the sense of increasing inferred accretion 



rate (e.g. Mendez & van der Klis 1999). These results 
are also consistent with the predictions of TLM. Recently, 
Cui (2000) proposed that the disappearance of kHz QPOs 
may be explained by the "disengagement" between the 
magnetosphere and the Keplerian disk. Campana (2000) 
suggested that the vanishing of the magnetosphere may led 
to the stopping of the kHz QPO activity. This could be 
also the natural result of TLM since the Keplerian oscilla- 
tor has no chance to enter the rotating frame of reference 
at all if it losts direct contact with the magnetosphere. 

We note that the similar QPO frequencies in the range 
of 1-400 Hz have been observed in X-ray millisecond pul- 
sar SAX J1808.4-3658 (Wijnands & van der Klis 1998a) 
and Galactic black hole candidates GRO J1655-40 (Remil- 
lard et al. 1999a), GRS 1915+105 (Morgan, Remillard & 
Greiner 1997), XTE .11550-564 (Remillard et al. 1999b) 
and XTE J1859+226 (Cui et al. 2000). It is still not clear 
whether the 67-400 Hz QPOs observed in these objects 
are identical to the kHz QPOs in neutron star LMXBs. 
But if we assume that these hundred-hertz QPOs are sim- 
ilar to the lower kHz QPOs, the correlations of their fre- 
quencies of low-frequency QPOs (l-20Hz) and break fre- 
quencies with the frequencies of these 'kHz QPOs' seem 
to be similar as found for atoll sources (see Psaltis et al. 
1999a; Wijnands & van der Klis 1999). In fact, if we plot 
the QPO frequency data for SAX J1808.4-3658 in Fig- 
ures 3 and 4, these points will locate in the lowerleft parts 
of both figures but are still consistent with the correla- 
tions for atoll sources. Their positions in these two figures 
also support that SAX J1808. 4-3658 is similar to the low- 
luminosity LMXBs with low accretion rate (Wijnands & 
van der Klis 1998a). Its estimated magnetic field strength, 
S < (2 - 6) X 10* Gauss (Wijnands & van der Klis 1998b), 
is indeed similar as that for atoll sources. For black hole 
candidates, a similar transition layer as in LMXBs may 
also exist between a geometrically thin clod disk and a ge- 
ometrically thick hot disk (or accretion-dominated accre- 
tion flow; see Narayan & Yi 1994) due to the sub-Keplerian 
rotating feature of the inner hot disk. Because of the lack 
of magnetosphere near the black hole, we can only predict 
three characteristic frequencies (ffe, Vy and Vbreak) from 
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TLM. This may help us to understand why we did not ob- 
served QPO frequencies larger than 400 Hz for black hole 
candidates and why their power density spectra and cor- 
relations between low and high QPO frequencies are more 
similar as those for atoll sources than for Z sources (van 
der Klis 1994, Psaltis et al. 1999a). A detailed compari- 
son of these features, however, is beyond the scope of this 
paper. 

We should also mention that as in other QPO models, 
TLM requires blob be lifted from the inner part of Ke- 
plerian disk. It is still unclear which mechanism could 
lead to such lifting. In addition, most theoretical predic- 
tions of TLM are made by assuming the angular velocity 
of magnetosphere SI is constant. This may not be the case 
in reality. Considering a non-constant SI will lead to the 



complexity of the dispersion relation and therefore some 
different characteristic oscillation frequencies. Moreover, 
the lack of explicit knowledge about the diffusion process 
in the transition layer may result in the uncertainty in es- 
timating the break frequency. In spite of these questions, 
TLM should be considered a very competitive model for 
QPOs in X-ray binaries. 

I am very grateful to Dimitrios Psaltis, Eric Ford, 
Michiel van der Klis, Rudy Wijnands, Mariano Mendez 
and Lev Titarchuk for kindly providing me the data files. I 
thank Wei Cui and the anonymous referee for helpful com- 
ments, and thank Lev Titarchuk, Wenfei Yu and Cheng- 
min Zhang for stimulating discussions. 
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